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Abstract 

We address the controversial hot question concerning the validity of the loose-coupling versus 
the lever-arm models in the actomyosin dynamics by re-interpreting and extending the washboard 
potential model proposed by some of us in a previous paper. In the new theory, a loose-coupling 
mechanism co-exists with the deterministic lever-arm model. The synergetic action of a random 
component, originating from the harnessed thermal energy, and of the power-stroke generated by 
the lever-arm classical mechanism is seen to yield an excellent fit of the set of data obtained in 
T. Yanagida's laboratory on the sliding of Myosin II heads on actin filaments under various load 
conditions. Our theoretical arguments are complemented by accurate numerical simulations, and the 
robustness of theory is tested via different combination of parameters and potential profiles. 
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1 Introduction 

In 1985 Yanagida et al. [J described a sophisticated experimental set-up to measure the distance traveled 
by an actin filament interacting with a Myosin II head 1 during a complete ATP cycle. Under low load 
conditions on the actin filament, they could observe traveled distances of up to 60 nm. This was in strong 
disagreement with other experimental evidences (see [2], PI and references therein), indicating instead 
distances ranging between 4 and 10 nanometers, such smaller range of traveled distances is coherent with 
the Lever-arm theory (see 0] and [H]) that is still widely believed to account for the generation of the 
force responsible for the actin filament sliding. The level-arm theory can be essentially schematized as 
follows (more details can be found in 6 , where a graphical animation is also present): 

- The binding of an ATP molecule on the catalytic site of myosin head produces the detachment of 
the head from the actin filament. 

*Work partially supported by MIUR (PRIN 2003) and by Campania Region. 
1 Hereafter we shall use the word "myosin" as an abbreviation of "Myosin II" . 
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- The ATP hydrolysis, that follows within about the next 10 ms, changes ATP to ADP.Pi with 
simultaneous release of energy in the amount of about 20 k B T. ATP hydrolysis also generates a 
sensible rotation of the neck of the myosin head. 

- The transition M. ATP — ► M. ADP. Pi is accompanied by an increase of affinity of this complex for 
actin, which enhances the chance of myosin to get binded on an actin site. Since such binding is 
reversible, the myosin head can visit more than a single actin site. 

- While myosin head is binded at an actin site, its neck may suddenly switch back to its original 
position (this is the so called "power stroke" ) with the successive release of the phosphoric radical 
Pi and the occurrence of conformational changes around the myosin coiled coil 2 . This process is 
reversible i.e. phosphoric radical Pi can re-establish the complex M. ADP. Pi. In the absence of 
such recombination, the above-mentioned conformational changes are very likely to generate the 
sliding of the actin filament with the release of ADP molecule after about 2 ms. 

- The myosin head stays then binded on the actin site till the arrival of a new ATP molecule, which 
starts the whole process afresh. 

Summing up, as far as reference is made exclusively to the mechanical phenomenon of the generation 
of the force responsible for the movement, the Lever-arm theory is strictly deterministic: Each ATP cycle 
generates one single power-stroke that causes a sliding of constant length and preassigned direction of 
the actin filament. Within such framework, a tight coupling between ATP cycles and protein movements 
is envisaged. 

A dramatically different conclusion appears to be required by the data of Yanagida and co-workers as 
indicated in the sequel. Indeed, aiming at a continuous efforts towards increasingly accurate measurements 
of the traveled distances during an ATP cycle, they were able to achieve impressive improvements of the 
technological set-up, including the design and construction of "home made" highly sophisticated devices. 
As described in [8], by exploiting such a technology, they proceeded as follows: A myosin head was 
attached to the tip of a glass microneedle and placed near an actin filament that had been previously 
immobilized on a microscopy slide by means of optical tweezers. The deflections of the needle with respect 
to its resting position were then measured and recorded. From the traces thus obtained the following 
three features emerged that are in evident disagreement with the tight coupling theory: 

1. The overall displacement traveled by the myosin head is not constant; it can be as large as 30 nm. 

2. This distance is the sum of a random number of single "steps", the amplitude of each of which 
equals the distance (~ 5.3 nm) between two successive actin monomers. During the time elapsing 
between two successive steps, the myosin head randomly jitters around an equilibrium position. 

3. Steps mainly occurs in a fixed "forward" direction, although some of them occasionally take place in 
the opposite "backward" direction. Hereafter, forward steps will be taken as positive and backward 
steps as negative. The overall displacement is thus the algebraic sum of the number of performed 
forward and backward steps. 

Such evidence contrasts the one-to-one relation hypothesized between ATP hydrolysis and the occurrence 
of the mechanical event consisting of the power stroke in the myosin head. Furthermore the observation 
of the existence of random elements leads one to conjecture that a significant role could be played in 
such context by the thermal agitation of the environmental molecules of the watery solution in which the 
involved proteins are embedded. This is the motivation for the assumption of the existence of a loose 
coupling between ATP cycle and actomyosin dynamics (See, for instance, [§])• 

While referring to |l()j for a lucid outline of the origin of the controversy existing among the supporters 
of the loose coupling mechanism and the community of those faithful to the Lever-arm theory, and 
therefore also the tight coupling vision, in the present paper we purpose to resolve such a controversy 

2 The nature of such conformational changes has not been yet fully understood. With reference to kinesin, a possible 
mechanism is described in where its validity also for the myosin family is conjectured. 
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by proposing a model that integrates the Lever-arm theory with certain components of random nature. 
Specifically, via the comparison of our theoretical results with those yielded by the experiments, we shall 
test the assumption that during the time interval elapsing between ATP hydrolysis and the final release of 
the phosphoric radical, the position of the myosin head is determined as the result of a dynamical process 
that includes a macroscopic non deterministic force responsible for a non-zero average net displacement. 
As a consequence, the coupling between the ATP cycle and the mechanical effects should appear to be less 
rigid, and the actomyosin dynamics should definitely exhibit variability features of the kind mentioned 
in the above^^El items. 

The conjectured random dynamics above outlined will be described in details in Section 2 with 
a specific reference to an earlier paper Here we limit ourselves to showing that, within such a 

framework, it is possible to handle a key point of the mentioned controversy. Indeed, in the essential 
relevance of the role played by the length of the myosin head is stressed, since the Lever-arm theory 
makes the sliding distance less than, and somewhat proportional to, such length. This would for instance 
imply that reducing the length of the myosin neck to a half should reduce to a half the sliding distance, 
which appears to be confirmed by the experiments in |12| . On the contrary, experiments performed by 
the Yanagida group only show slight changes of the overall displacement of the myosin head even after 
complete removal of its neck |10) . 

In order to reconcile these evident discrepancies, we shall assume that the displacement X of the 
myosin head during the ATP cycle can be represented as follows: 

X = rX R + dX D (1) 

where Xr denotes the displacement induced by the random force and Xp the displacement generated by 
the power stroke, and where r and d are constants, each equal to or 1. Note that the case r = and 
d = 1 yields the Lever-arm theory; instead, setting r = 1 and d — depicts a purely random situation, 
in the absence of any sliding due to the power stroke. Finally, the case r = d = 1 leads to an integration 
of the two theories. With the choice r = d — 1, the controversial results related to the length of the neck 
of myosin head can be overcome by assuming that the random displacement is a few times larger than 
the deterministic one. This assumption implies that only slight changes of the total distance traveled 
by myosin head would be observable when the contribution Xjj due the power stroke is made somewhat 
smaller by shortening the length of the myosin neck. 

2 The model 

Let L denote the distance between each pair of neighboring actin monomers. As suggested in recent 
literature (see |13| and [14| ) in our computations we shall take L = 5.5 nm. In addition, we shall assume 
that the magnitude of the sliding induced by the power stroke equals that of a step during the random 
phase, and thus take Xd — L. 

Our conjectured relation Q with r = d = 1 is preliminarily supported by the bar chart in Fig. 4c) 
of 8 showing that at least one step is performed by each and every myosin head during the rising phase 3 . 
The second column of Table ^ shows the heights of the columns of the mentioned bar chart, whereas the 
first column indicates net step numbers, i.e. the integral part of the ratios X/L. The quoted experiment 
was performed under low load conditions by means of microneedles having stiffness less than 0.1 pN/nm. 

As shown by the first two columns of Table ^ one sees f° r instance that out of 66 observed myosin 
heads (all attached to a glass microneedle of stiffness less than 0.1 pN/nm) 14 performed a unit net steps, 
21 a net step number equal to 2, etc., to conclude that 3 of them have performed a net step number 
equal to 5 implying total displacement of about 30 nm. Columns 3 and 4 of Table ^ show that the net 
step number \Xr/ ' L\ of the random component is well fitted by a Poisson distribution with parameter 
fj = 1.5 given by the ratio of the total number of performed net steps (99) to the number of considered 
myosin heads (66). 

3 The rising phase is the time interval starting with the hydrolysis of the ATP molecule and ending with the final release 
of the phosphoric radical Pi. 
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Table 1: Observed distribution of net number of steps performed by myosin 
heads as indicated in [S]- Here [_^7-^J is the total step number during the entire 
rising phase, whereas \_Xr/L\ = [X/L\ — 1 denotes the total ( algebraic) step 
number during the random part of the rising phase. 



[X/L\ 


Observed frequency 


[Xr/l\ 


Theoretical frequency 


1 


14 





15 


2 


21 


1 


22 


3 


18 


2 


17 


4 


10 


3 


8 


5 


3 


4 


3 


6 





5 


1 


total 


66 


total 


66 



The agreement between experimentally observed frequencies and those predicted via the Poisson dis- 
tribution, jointly with the rarity of the backward steps, leads one to conclude that the "dwell time" 
(namely the time interval elapsing between two successive steps) should be, to a good approximation, ex- 
ponentially distributed. Such conclusion is experimentally supported by the data summarized in Fig. 4b) 
of 8 leading to the estimate of about 5 ms for the mean dwell time. It must be pointed out that in the 
mentioned estimation the dwell time of all performed steps have been included, whereas our conjecture 
on the exponential distribution of dwell time only refers to the steps performed during the random phase. 
However, such a diversity does not lead to inconsistencies if one assumes that the dwell time between the 
last step due to the random phase and the step due to the deterministic power stroke is exponentially 
distributed as well. 

Our facing sequences of rare events with exponentially distributed interarrival times is strongly sug- 
gestive of a first-exit problem out of an interval for continuous Markov processes possessing an equilibrium 
point sufficiently far from at least one of the end points of the diffusion interval (See, for instance, 

On the ground of all foregoing considerations, with reference to the random rising phase we are lead 
to construct a model for the actomyosin dynamics that is based on the following assumptions: 

(i) The complex M.ADP.Pi + Energy is viewed as a point-size particle moving along an axis X on 
which the abscissa x denotes the displacement of the particle from the starting position. The 
positive direction of X is that of the forward steps of myosin head. 

(ii) The particle is embedded in a fluid. Hence, it is not only subject to a dissipative viscous force 
characterized by a drag coefficient /?, but also to microscopic forces originating from the thermal 
motion of the fluid molecules. On account of the fluctuation-dissipation theorem, such microscopic 
forces can be macroscopically described by means of a Gaussian white noise having intensity 2f3k B T, 
where k B is Boltzmann constant and T the absolute temperature. 

(iii) The global interaction of the particle with the actin filament is synthesized in a conservative unique 
force deriving from a potential U(x). The structure of the actin filament suggests that U(x) be a 
periodic function with period equal to the distance L between pairs of consecutive actin monomers: 

U(x) = U{x + rL), VreZ. (2) 

Henceforth we shall denote by La (0 < La < L) the minimum of U(x), assume U(La) = and 
denote by Uq :— U(0) = U(L) the depth of the potential well, namely the maximum of U(x). 

(iv) The particle's dynamics is described by Newton's equation in which the total acting force is the 
sum of two terms: the first term is a deterministic force generated by the potential U(x) and by 
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the viscous force, while the second term is the random force due to the presence of the Gaussian 
white noise. 

(v) One more (constant) force F{ > acts on the particle. Here we assume that this force is generated by 
a process that finds its origin in a part of the energy possessed by the particle, and take FjL <C Uq. 

Summing up, we are assuming that the complex M.ADP.Pi + Energy can be looked at as a Brownian 
particle subject to a tilted potential V(x): 

V(x) = U{x) - Fx (3) 

where U(x) has been indicated in J5J), F :— Fi ~ F e , with F e > 0, is a constant external force, possibly 
applied from the outside by the experimenter. In the experimental conditions the height of the potential 
wells is Uo < 100 pN-nm, the period of the potential is L = 5.5 nm, the particle mass is rn = 2.2 • 10~ 22 kg, 
the drag coefficient is j3 — 90 pN-ns/nm and the environmental temperature is T — 293 K. Therefore, 
Reynolds's number is much less than 1 (see, also, ^Sj), so that the inertial term of the equation of motion 
can be disregarded. In conclusion, the overdamped equation describing the movement of the particle is 
the following Langevin equation: 

2knT m (4) 

where A(t) is a zero-mean white Gaussian noise with unit intensity, ( ' ) denotes time derivative and ( ' ) 
space derivative. 

Within this framework, this idealized Brownian particle randomly moves around an equilibrium point 
located at a minimum La of U (x) . Whenever it exits the current potential well, we conventionally 
say that the corresponding myosin head has made a step, in the forward or in the backward direction 
according to where the exit has taken place. Hence, we are facing a first-exit problem of the Brownian 
particle from the endpoint of the current potential well. Taking into account the above assumptions, we 
conclude that the distribution of the first-exit time of the Brownian particle from the potential well is 
exponential, since: 

- the process described by equations J5J, J3J) and Q possesses an equilibrium point at the minimum 
La of U(x)- 1 

- the time for the particle to travel the distance La in the presence of the only force due to U(x) is 
(see for instance, ^3) r = @L A /Uo; 

- the standard deviation of the Gaussian steady-state distribution of the process modeling the parti- 
cle's motion is a = y/ (k B T/ (3) ■ r/2 = La/V^uq, where we have set uq — U~o/k B T; 

- the ratio I a '■= La/<t = V^uo falls well within the interval (2,4-^/2) for all choices of ito G [2, 16] to 
which we shall refer in the foregoing, so that the exponential approximation for the first-exit time 
is valid [15] . 

In the sequel, we shall exploit the known formulas ^Hl f° r the probability p that the particle exits 
from the current potential well to reach next well, 

P = l + exp(-FL/k B T)' ^ 
and for the mean first-exit time fj, from a potential well: 
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3 Determination of F{ and Uq 

Hereafter, we shall view as constants the parameters /3 and T that characterize the thermal bath and 
the period L of the potential U(x), their values being specified as in Section [21 Hence, the quantitative 
specification of the model described by Eqs. (J2J, J3J) and (0J requires that numerical values be attributed 
to three more parameters: the depth Uq of the potential well, the position La of the minimum of U(x) 
in (0, L) and the internal force Fi. The numerical specification of these parameters can be performed 
after the function U(x) has been chosen. We shall preliminarily take U(x) as a symmetric {La = L/2) 
saw-tooth potential, as defined in Table El Successively, we shall test the robustness of our model by 
assuming alternative potential functions, henceforth called "potential profiles". It should be noted that 
Fi is somewhat related to the largest force that myosin is able to endogenously generate. Here, we shall 
not attempt to provide any biological justification of its genesis. We limit ourselves to pointing out that 
elsewhere JH1, where similar experiments were performed on Myosin VI, it is conjectured that the week 
binding between actin and myosin is a source of distortion of the geometry of the two helixes in the actin 
filament. Such a distortion exposes the hydrophobic region of actin to myosin head, thus generating a 
tilt of the potential, and hence a constant force Fi. The existence of a tilt of the potential, conjectured 
in has been successively supported by means of simulations in |20| where it is shown that in the 
absence of such a tilt experimental available evidence on the myosin motion in the presence of contrasting 
applied loads cannot be accounted for by any of the other models therein considered. Within our strictly 
phcnomcnological framework the existence of this force Fi is supported by Eq. JSJ showing that in the 
absence of external applied forces (i.e. F e = 0) steps on either directions would be equally likely unless 
Fi =/= 0, in contrast with the experimental evidence on the high degree of directionality exhibited by 
motion of the myosin head. 

To proceed along the quantitative specification of our model in a way to be able to attempt the fitting 
of available experimental data, the values of Fi and of Uq must be specified. This will be done by making 
use of the available experimental data (2H shown in Tableland in Table Of. From them it is evident that 



Table 2: For three conditions of the applied load C the 

recorded numbers hf of forward steps, h b of backward steps, Table 3: Recorded dwell times fi 

the total number of steps and the percentage p of forward for different values of the applied 

steps are listed. load C. 



C (pN) hf h b hf+h b p—- — 

rif + n b 

] 0.0, 0.5] 54 9 63 0.8571 



] 0.5, 1.0] 40 9 49 0.8163 



] 1.0, 2.0] 29 19 48 0.6042 



the frequency p of forward steps decreases as the applied load increases, while the dwell times increase 
with the load. Since the external space-dependent force accounting for the applied load is not included 
in our above formulated model, one more assumption is necessary: 

(vi) The external constant force F e acting on the myosin head equals the value C obtained by the 

4 Noto that loads and dwell times in TablelBlmust be viewed as averages to which confidence intervals are associated. For 
instance, the load 0.046 pN is the result of all measurements for which the product of the stiffness of the glass microneedle 
times the distance traveled by the myosin head falls around 0.046 pN. The corresponding dwell time 5.3 ms is to be viewed 
as the arithmetic average of the dwell times recorded during these measurements. 



C (pN) 


pu (ms) 


0.046 


5.3 


0.190 


5.7 


0.300 


6.0 


0.470 


7.1 


0.690 


8.9 


0.830 


6.2 


1.240 


11.1 


1.890 


11.0 
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clastic force C (x) (given by the product of the microneedle stiffness times the distance traveled by 
the myosin head) at the end of the rising phase and its direction is opposite to that in which motion 
occurs. 

In other words, in our model the presence of applied loads in the experimental setup is expressed by 
means of the variable F e . 

For some fixed values of internal force Fi , use of Equation JSJ) has been made to calculate the theoretical 
probabilities of the particle's exit from the current potential well to the next well (i.e. the analogue of 
the forward step probabilities) as function of F e . The results are shown in Figure ^ where eight realistic 
values of Fi have been chosen in the interval [l.OO pN, 1.90 pN] . Vertical lines indicate the 3 load 
intervals of Table |3 whereas horizontal lines indicate the 3 corresponding recorded frequencies. We see 
that for internal forces Fi — 1.00 pN and Fi = 1.90 pN the plotted curves do not meet the requirement of 
leading to the experimentally recorded frequency p = 0.8571, whatever value F e is chosen within interval 
[0 pN, 0.50 pN] . Similarly, for Fi = 1.55 pN no value of the computed probability equals the frequency 
p = 0.8163 for F e ranging in [0.50 pN, 1.00 pN] . Instead, all remaining 5 curves referring to values of Fi 
ranging from 1.60 pN to 1.80 pN in steps of magnitude 0.05 pN, from below upward, are in agreement 
with the experimental values of Table [21 Hence, the interval of values for Fi to be selected in order to 
secure the fitting of the experimental data is [1.60 pN, 1.80 pN] . 

1.0 i 1 




0.3 - 



0.2 -I 1 1 1 ^ 

0.0 0.5 1.0 1.5 2.0 

External applied force (pN) 

Figure 1: The probability p that the particle exits from the current potential well to 
reach next potential well is plotted as function of the external applied force F e for 
various values of internal force Fi (pN). Temperature and potential period are taken 
as T = 293 K and L = 5.5 nm, respectively. 

We now come to the estimation of the last parameter, Uo, i.e. of the depth of the potential well. This is 
done by exploiting Eq.[B]in which the left-hand side is viewed as the function /i = [i(Fi, F e , Uq; L, /?, T) and 
fixed to 5.3 ms, that (see TableEJ corresponds to the smallest applied load 0.046 pN 5 . Since temperature 
T, period L of the potential V(x) = U(x) — (Fi — F e )x, and drag coefficient j3 are specified, if we take 
F e = 0.046 pN Eq. @ makes Uo an implicit function of Fi. 

5 We choose the smallest load to minimize the error due to the approximation of elastic force C(x) by F e . 
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Note that as Fi increases (i.e. as the potential tilt increases) the mean first-exit time \i decreases. 
Hence, in order to keep \i constantly equal to 5.3 ns, while Fi ranges in the interval [l.60 pN, 1.80 pN] , 
depth Uq must be taken as a monotonically increasing function of Fj. From /z(1.60, 0.046, Uo) = 5.3 and 
/x(1.80, 0.046, U ) = 5.3, Eq. © yields U « 15.632 k B T and U « 15.755 fc B T, respectively. 

In order to test the agreement of our model with the experimental values of mean dwell times for 
the various loads (see Table OJl , we make use of Eq. © to determine \i as a function of F e for the pairs 
(1.60 pN,15.632 k B T) and (1.80 pN,15.755 k B T), involving the extrema of the determined values for Fi 
and Uq. The obtained values are showing Figure where the corresponding experimental mean dwell 
times are also indicated. Obviously, all other pairs of admissible values of Fi and Uo lead to curves lying 
inside of the above two pairs. Inspection of Figure |2 jointly with the magnitudes of the confidence 



12.0 




0.0 0.2 



0.4 0.6 0.8 1.0 1.2 1.4 
External applied force (pN) 



1.6 1.8 2.0 



Figure 2: First-exit time is plotted as a function of the external applied force 
F e . The chosen values of Fi (pN) and Uo (k B T) are indicated for each curve. Dots 
represent the experimental dwell times of Table Drag coefficient has been taken 
as 90 pNns/nm, whereas temperature and period of the potential are the same as in 
Figure d 

intervals associated to each experimental value ^1], shows that for small values of F e the agreement 
between experimental data and theoretical predictions is good. Indeed, up to the first 3 values of Tabled 
the agreement is excellent, to become more than acceptable for larger values of F e up to 0.83 pN. The 
discrepancies shown by the remaining two experimental data are a consequence of the crude assumption 
(vi) above by which the applied load is assumed to be strictly parallel to the direction of motion. This is 
acceptable for reasonably small loads, but certainly unrealistic for large loads: In this case a significant 
orthogonal component should be expected, whose effect is alike to an increase of the depth of the potential 
well Uo- Nevertheless, in the interval between the last two load values of Table a qualitative similar 
behavior of theoretical curves and experimental values is present. 
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4 Net step number distribution 

In the present Section we shall implement the model described by Eqs.J2J l@J in order to obtain, via a 
suitable simulation procedure, the distribution of the net number of steps performed by the myosin head 
during the time interval elapsing between the instant when the ATP molecule is hydrolyzed and the final 
release of the phosphoric radical, i.e. during the random part of the rising phase. The results of our 
simulations will then be discussed with reference to the experimental data of Table ^ We recall that 
such data refer to 66 myosin heads that have globally performed 99 net steps. 

Our simulation procedure is based on a discretized version of Eq. and on a routine for generating 
Gaussian pseudorandom numbers in a way to determine step by step the positions achieved by 66 Brow- 
nian particles all originating at xq at time 0. The traveled distances at the end of the individual random 
rising phases are recorded, which finally leads one to the determination of the 66 net step numbers. Such 
procedure has been implemented on an IBM SP4 parallel supercomputer and repeated 1600 times to 
obtain a reliable statistics. 

In order to solve Eq. Q numerically, one has to specify preliminarily initial position xq of the Brownian 
particle and duration O of each sample path. We shall safely take xo = L/2 since, whatever the actual 
initial position of the particle inside of the potential well, the its relaxation time is much smaller than the 
mean first-exit time. The specification of sample path duration is somewhat more involved because is 
a random variable of unknown distribution. Nevertheless, we can appeal to the approximate exponential 
distribution of the first-exit time, motivated by the depth of the potential well, to assume that O is 
approximately gamma-distributed, with probability density v) where /z is the mean first-exit time 
from a potential well and v is the mean total exits of the Brownian particle. While \x is obtained via 
Eq. 10, our estimate v of v is obtained by using the data of Tableland Eq. J5j|: 

- - 99 

V ~ 66(2p- 1)' 

Finally, for each specified potential well Uq, the size At of the time parsing in Eq. JJJJ is determined by 
progressively reducing it until the obtained distribution becomes appreciably invariant. For an immediate 
comparison, Table 0] shows the distribution obtained via simulation for a potential Uq = 5 k B T and a 
parsing time At — 0.25 ns jointly with the experimental distribution. The agreement between observed 
and numerical frequencies is more than satisfactory. 

The chosen value of Uo is motivated by a twofold consideration. First of all, it is large enough to 
imply that the first-exit time distribution is approximately exponential. This is indeed supported by the 
numerically evaluated variation coefficient that has been found to be 0.96. Second, and most relevant 
consideration, is that it is reasonable to conceive that, under the assumed overdamped regime, the net 
step number is insensitive to the depth of the potential well. Indeed, the forward exit probability p 
given by @, under the fixed environmental temperature, only depends on the energy FL, namely on the 
difference of potential V(x) over one period, thus being independent of the depth Uo- Table evidently 
supports such conclusion. Indeed, it indicates that, for instance, by doubling the depth of the potential 
well Uo, the net step frequencies are not affected significantly, even for different choices of internal forces 6 . 

It should be explicitly remarked that the parsing parameter At must be determined with specific 
reference to the magnitude of Uq- Indeed At must be much smaller than the relaxation time of the force 
— U'(x) that is proportional to f3L 2 /Uo. Furthermore, it must be such as to cope with the random forces 
due to the thermal bath. Therefore, the mean square displacement per unit time should remain constant 
as the depth Uq of U(x) is made to change, which implies an inverse dependence of the magnitude of At 
on the square of magnitude of the potential well. This is the motivation for the indicated choices of At 
in Table El 

6 The presence of non-integer numbers in Tabic l4| Island Islis a consequence of our adopted estimation procedure. The 
half-width of the related 95%-confidence intervals has been seen never to exceed 0.2. For comparison reasons we have not 
rounded out the raw numbers to the nearest integers. 
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Table 4: Net step frequency distributions for Uq — 5 k B T and for three values of Fi chosen 
in the admissible interval. Here F e — 0, xo — L/2 and At = 0.05 ns. Other parameters have 
been chosen as in Figure |21 The last 5 rows show the coefficient of variation, probability of 
a step forward and mean exit time, all obtained via simulations. The theoretical values of p 
and //, obtained via Eqs. (0) and JJjJ, are also indicated. 



[X/L\ 


Observed frequency 


Numerical frequency 






Fi = 1.70 pN 


Fi = 1.75 pN 


Fi = 1.80 pN 


-1 





2.05 


1.90 


1.76 





14 


13.79 


13.63 


13.88 


1 


21 


20.72 


20.83 


20.93 


2 


18 


16.41 


16.43 


16.47 


3 


10 


8.53 


8.59 


8.54 


4 


3 


3.14 


3.26 


3.16 


5 





0.93 


0.93 


0.88 


CV 


0.966 


0.964 


0.963 


P 


0.909 


0.914 


0.920 


P 


0.910 


0.915 


0.920 


~P (ns) 


1235.057 


1207.288 


1180.852 


fi (ns) 


1233.784 


1206.048 


1178.695 



5 The role of potential forms and asymmetries 

In |llj a preliminary version of the model considered here was discussed with reference of a parabolic 
potential, henceforth called "parabolic profile". The evaluation of the robustness of the model was, 
however, postponed to a successive investigation. This task is the object of the present Section in which, 
in addition to the case of saw-toot profile, three more profiles will be considered: parabolic, cosine- like and 
Lindner- type. (See Table EJ- A plot of the considered potentials over one period are shown in Figure |3| 
whereas Figure 01 refers to the corresponding generated forces. Lindner-type potential has been indicated 
for two values of parameter 5. It is not difficult to see that S — > yields the cosine profile, whereas the 
potential flattens down in the middle as 6 increases. 

The independence of the exit probability p of the potential's profile implies that the net step distri- 



Table 5: Net step frequency distribution. Here F e — and xq — L/2. Other parameters are chosen as 
in Figure 13 Parsing steps are chosen as follows: At = 0.1 ns for Uq = 4 k B T, and At — 0.025 ns for 
Uq = 8 k B T. The indicated values of Fi belong to the admissible interval. 





Fi = 1.70 pN 


Fi = 1.75 pN 


Fi = 1.80 pN 


[Xa/L\ 


U = 4 k B T 


U = 8 k B T 


U = 4 k B T 


U = 8 k B T 


U = 4 k B T 


U = 8 k B T 


-1 


2.11 


2.05 


1.94 


1.90 


1.85 


1.81 





13.56 


13.96 


13.54 


13.76 


13.48 


14.00 


1 


21.13 


20.26 


20.92 


20.26 


21.30 


20.46 


2 


16.70 


16.11 


16.93 


16.29 


16.98 


16.09 


3 


8.43 


8.58 


8.53 


8.61 


8.40 


8.56 


4 


2.95 


3.40 


3.01 


3.53 


2.91 


3.49 


5 


0.75 


1.12 


0.78 


1.15 


0.76 


1.11 
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Table 6: Potentials' profiles. 



Saw-toot 



U s (x) 



-Up 
La 
Uo 



L-L 



(x — La), < x < La 
(x - L A ), L a <x<L 



Parabolic 



Up{x) = jxji V i 



Cosine 



Uo(x) = f 



-X ] J 



Lindner 



U L {x) = 



,26 



- 1 



cos I —x ) + 1 



- 1 



Table 7: For each potential pro- 
file the depth Uq of the poten- 
tial well is indicated. Here Fi = 
1.75 pN, F e = 0.046 pN and the 
other parameters are the same as 
in Figure [3 



Type 


Uo (k B T) 


Saw-toot 


15.723 


Parabolic 


15.043 


Cosine 


13.944 


Lindner (8 = 0.1) 


13.942 


Lindner (8 = 0.5) 


13.918 


Lindner ((5 = 1) 


13.851 


Lindner (8 = 2) 


13.697 


Lindner (8 — 3) 


13.611 


Lindner (8 — 4) 


13.591 


Lindner (8 = 5) 


13.604 


Lindner (8 = 10) 


13.749 



bution is profile- independent as well. This is also evident from Table |H1 where the net step distributions 
are reported for each of the above-considered four potential profiles. 

Next task is to pinpoint the effects of the potential profiles on the mean first-exit time. To this 
purpose, we refer to Table [3] showing that the mean dwell time in lowest load condition is \i = 5.3 ms. 
After choosing Fi = 1.75 pN, we then make use of Eq.JU) for each and every one of the four considered 
potentials imposing that the left hand side equals such value of /i. By iterated numerical integrations, 
for each potential profile the corresponding value of Uo is finally determined. The result are listed in 
Tabled where the effect of the parameter 8 in Lindner- type profile has been detailed. 

The behavior of mean first-exit time n as a function of the magnitude of the external force is shown 
in Figure [3] All considered cases of Table [7| lead to graphs falling in the region bounded by the lowest 
and the highest curves. We point out that changing the potential profiles never yields mean first-exit 
time changes exceeding ten percent. Hence, we are led to conclude that the depth Uo of the potential 
well can be tuned to acceptable biological values by a suitable selection of the potential profile, without 
affecting appreciably the value of the mean first-exit time. For instance, switching from saw-toot to 
Lindner-type potential with 8 = 2, lowers Uo by more than 2 k B T. (See Table 0). It is thus conceivable 
that a variety of potential profiles exists such that the height of the potential well can be further lowered, 
in a way to switch from about 15 k B T as indicated in [20] to about 5 k B T as suggested in jEj. The 
quantitative analysis and comparison with available data has been performed under the assumption of 
rigorous symmetries exhibited by the profiles of the potentials generating the periodic force acting on 
the Brownian particles. It is, however, presumable that the complex biological reality underlying the 
observed motion of the myosin heads may require to relax the rigorous symmetry assumption. To test 
the effect of symmetry breaking, we take into consideration the saw-toot potential Us(x) of Table Eland 
make it asymmetric by taking La L/2, namely by shifting in either direction the point of minimum of 
the potential. Thus doing, the introduced asymmetry should affect the motion of the Brownian particles. 
While the probabilities of exit from the current potential well are insensitive to the potential's profile, 
and hence also to its asymmetry, the mean first-exit time [i is clearly affected by it, as shown by Eq. ©. 
The quantitative dependence of fi on the potential's asymmetry is indicated in Figure where on the 
abscissa the external force F e is indicated. Each curve is characterized by two parameters: the point La 
of minimum and the depth Uo of the potential. From top to bottom, the first curve is characterized by 
asymmetry La = 2.0 nm, implying a somewhat steeper rise of the potential in the backward direction. 
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Figure 3: Plots of the potentials as function 
of the position within each pair of consecutive 
monomers. 



Figure 4: Conservative forces originated by the 
potentials of Figure 



Such asymmetry is mitigated in the next curve (La = 2.5 nm). The third curve, indicated in bold, is 
used as a comparison tool, as it refer to symmetric profile (La = L/2). The remaining two curves are 
labeled by La = 3.0 nm and La = 3.5 nm, thus representing mirror-image situation with respect to the 
zero-asymmetry case. Now the forward edges of the potentials are steeper. For each curve, the indicated 
value of Uo has been determined by imposing that /i = 5.3 ms when F e = 0.046 pN (lowest load), so that 
all curves originate at the same point. As Figure El shows, changes of the potential's asymmetry of 30% 
in either direction do not affect greatly the magnitude of the mean first-exit time. Finally, we remark 
that suitable choices of backward asymmetry (such La — 2.0 nm in Figure El improve the fitting of the 
experimental data even for larger load values. 



6 Concluding remarks 

The model of actomyosin dynamics discussed in the foregoing rests on the assumption that the total 
energy made available to the myosin head by the ATP molecule hydrolysis and by the thermal bath has a 
two-fold overall role: To produce the power stroke predicted by the lever-arm model and also to generate 
the kind of sliding of myosin head on the actin filament in the "loose coupling" mechanism originally 
hypothesized in [22 and then experimentally demonstrated in [5]. This is expressed mathematically via 
the representation of the displacement of a particle consisting of a combination of a deterministic part 
and of a random component. The latter is generated by the simultaneous presence of a washboard- type 
potential and a random force arising from thermal fluctuations. Our model has then been tested by 
making use of a set of data on the dwell times and step frequencies of myosin heads under various load 
conditions. We have shown that by a suitable tuning of the internal force and depth of the potential 
well, the theoretically calculated probability p and mean first-exit time /i of the representative Brownian 
particle are in good agreement with their biological counterpart. A second set of experimental data 
concerning the net step number distribution of myosin heads under low load conditions has then been 
exploited to show that the washboard potential used by us is able to reproduce such distribution within 
the mathematical analogy. To the rather small size of the experimental sample should be ascribed the 
3% discrepancy represented by the two net backward displacements predicted by our model. 

Next, the robustness of our model has been tested by inserting 4 different potential profiles in the 
Langevin equation of motion. The consequently performed calculations have shown that the mean exit- 
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Table 8: The first two columns list the net number of steps performed by myosin heads 
and the corresponding observed frequencies. The remaining four columns show the net step 
distributions, obtained via numerical simulations using the indicated potential profiles all 
possessing J7 = 5 k B T. In all cases F{ = 1.75 pN, F e = 0, xq = L/2, and At = 0.05 ns. 
Other parameter have been chosen as in Figure|3 Variation coefficients and mean-exit times 
/i, obtained via Eq. ©, are listed as well. 



[X/L\ 


Observed frequency 


Potentials' profiles 






Saw-toot 


Parabolic 


Cosine 


Lindner (5 = 2) 


-1 





1.89 


1.92 


1.90 


1.99 





14 


13.61 


13.69 


13.73 


13.81 


1 


21 


20.67 


20.69 


20.59 


20.46 


2 


18 


16.47 


16.42 


16.17 


16.20 


3 


10 


8.63 


8.60 


8.70 


8.64 


4 


3 


3.33 


3.29 


3.38 


3.34 


5 





0.98 


0.95 


1.05 


1.07 


CV 


0.964 


0.970 


0.984 


0.982 


M (ns) 


1206.048 


1403.504 


2194.007 


2270.677 



time depends on the external force in a way that is essentially insensitive to the chosen potential's profile. 
The chosen profile, instead, has been seen to play an essential role in that it critically relates the depth 
of the potential well to magnitude of the mean exit-time. A finer tuning of the mean exit-time is finally 
achieved by regulating the level of asymmetry of the potential's profile. 
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Saw-toot Parabolic Cosine Lindner- - - ■ 

9.5 n 




External applied force (pN) 

Figure 5: For each potential profile the mean first-exit time fi is plotted as a function 
of the external force F e . The corresponding values of Uq, for each potential, are listed 
in Tabled We have taken Ft = 1.75 pN, while the other parameters have been chosen 
as m Figure 
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